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Abstract: In this paper we present a stability proof of model predictive control without stabi- 
lizing terminal constraints of cost which are subject to unknown but measurable disturbances. 
To this end, a relaxed Lyapunov argument on the nominal system and Lipschitz conditions on 
the open loop change of the optimal value function and the stage costs are employed. Based 
on the special case of sensitivity analysis, we show that Lipschitz assumptions are satisfied if a 
sensitivity update can be performed along the closed loop solution. To illustrate our approach 
we present a halfcar example and show performance improvement of the updated solution with 
respect to comfort and handling properties. 

Keywords: suboptimal control; predictive control; stability analysis; Lyapunov stability; state 
feedback; stabilizing feedback 
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1. INTRODUCTION 

Due to its simple structure, model predictive control 
(MPC) has become a well-established method for (sub)op- 
timal control of linear and nonlinear systems, see, e.g., 
Camacho and Bordons [2004] and Rawlings and Mayne 
[2009], Grime and Pannek [2011]. By means of this method 
an approximated closed-loop solution of an infinite hori- 
zon optimal control problem is computed in the following 
way: in each sampling interval, based on a measurement of 
the current state, a finite horizon optimal control problem 
is solved and the first element (or sometimes also more) of 
the resulting optimal control sequence is applied as input 
for the next sampling interval(s). This procedure is then 
repeated iteratively. 

Unfortunately, stability and optimality may be lost due to 
the trunctation of the infinite horizon. In order to ensure 
stability of the resulting closed loop, one usually imposes 
terminal point constraints as shown in Keerthi and Gilbert 
[1988] and Alamir [2006] or Lyapunov type terminal costs 
and terminal regions, see Chen and Allgower [1998] and 
Mayne et al. [2000]. A third approach uses a relaxed Lya- 
punov condition presented in Grime and Rantzer [2008] 
which can be shown to hold if the system is controllable in 
terms of the stage costs, cf. Grime [2009]. Additionally, 
this method allows for computing an estimate on the 
degree of suboptimality with respect to the infinite horizon 
controller, see also Shamma and Xiong [1997] and Nevistic 
and Primbs [1997] for earlier works on this topic. 

Here, we follow the third approach but extend its appli- 
cability for the case of parametric control systems and 
subsequent disturbance rejection updates. In particular, 
we impose an abstract update of the control law with 
respect to the measured disturbance or parameter. Then, 
if Lipschitz conditions on the open loop change on both 
the optimal value and stage cost function hold for such an 



update, we can utilize the relaxed Lyapunov condition for 
the nominal open loop solutions to show stability of the 
resulting closed loop. 

Thereafter, we focus on the update law for optimal controls 
via sensitivities, see also Zavala and Biegler [2009]. Using 
results from sensitivity analysis, we can show that if a 
sensitivity update can be performed along the closed loop, 
then the Lipschitz assumptions hold and stability of the 
disturbed system can be guaranteed. Note that due to 
abstracting from the form of the update, the proposed 
stability proof is not limited to sensitivities but may also 
applicable for other update methods such as realtime 
iterations or hierarchical MPC, see Diehl et al. [2005] and 
Bock et al. [2007] respectively. 

The paper is organized as follows. In Section 2 the problem 
formulation and the concept of practical stability are 
defined. In the subsequent Section 3 the stability proof of 
MPC subject to disturbances without stabilizing terminal 
constraints and costs is given. In Section 4, the additional 
assumptions required within the stability proof are shown 
to be satisfied in case of sensitivity based updates of the 
control law. Moreover, we present simulation results for 
a halfcar subject to disturbed measurements of both the 
state and the sensor inputs in Section 5. To conclude our 
paper, we draw some conclusions and give an outlook on 
future research. 

2. SETUP AND PRELIMINARIES 

The types of control systems we consider within this work 
are given by the dynamics 

x(n + 1) = f(x(n),u(n), w(n)) (1) 

where x denotes the state of the system, u the external 
control and w an disturbance which can be measured. 
These variables are elements of respective metric spaces 
(X, dx), (U,du) and (W, dw) which represent the state, 



control and disturbance space. Therefore, our results are 
also applicable to discrete time dynamics induced by a 
sampled finite or infinite dimensional system. For ease of 
notation, we introduce the abbreviation = di(x,y) 

for di € {dx,du,dw}- Both, state and control, are con- 
strained to be elements of subsets X C X and U C U . 
We denote the undisturbed or nominal state trajectory 
corresponding to an initial state x a e X, a control se- 
quence u = {u(k)) keX and a nominal disturbance sequence 
w = (w(k)) kex , 2 := {0, 1, . . . , N - 1} with N e N or 
I := No by x(-) = x u . w {-; xo). Similarly, the disturbed state 
trajectory which is subject to a disturbed initial value Xo, 
a control sequence u = (u(k)) keX the disturbance sequence 
w = {w(k)) keX is denoted by x(-) = 5%,uj("i#o)- Since in 
the presence of constraints not all control sequences u are 
admissible, we introduce U x (xq, w) as the set of all admis- 
sible control sequences u — {u(k)) keX which for a fixed 
disturbance sequence w satisfy f(x(k),u(k),w(k)) e X 
and u(k) e U for k e T. 

Our aim in this work is to stabilize system (1) at a 
controlled nominal equilibrium, that is a point x* such 
that there exists a control u* and a nominal disturbance 
w* satisfying f(x*,u*,w*) — x* . To this end, we consider 
a two stage feedback design. In the first stage we want 
to compute a static state feedback law for a given nom- 
inal disturbance sequence w — (w(n)) neNo which mini- 
mizes the infinite horizon cost functional J 00 (x,u,w) = 
~Y^ = Q^{x{n),u{n),w{n)). In this context, the stage costs 
I : X x U x W — > Kg arc continuous and satisfy 
£(x*,u*,w*) — and £(x,u,w) > for all u € U for each 
x 7^ x* and each «i ^ w*. In order to avoid solving the 
discrete time equivalent to a Hamilton- Jacobi-Bellman 
equation which is computationally intractable in most 
cases, we use a model predictive control (MPC) approach 
to approximate the desired solution. The resulting cost 
functional is given by 

N 

Jn(x, u, w) := ^ £(x(k; x),u(k),w(k)) (2) 

k=0 

where N e N>2 denotes the length of the prediction 
horizon, i.e. the prediction horizon is truncated and, 
thus, finite. Consequently, the obtained control sequence 
u*(-,x,w) itself is also finite. 

In order to retrieve an infinite sequence, one usually im- 
plements only the first part of this sequence /j,n(x,w) := 
u*(0,x,w), then the prediction horizon is shifted forward 
in time and the procedure is iterated ad infimum, cf., e.g., 
Camacho and Bordons [2004], Alamir [2006] or Rawlings 
and Mayne [2009]. 

Remark 2.1. Note that even the evaluation of both 
and J/v require knowledge of a sequence of nominal future 
disturbances. Such values may be obtained by forward 
measurements, e.g. cameras in a car or thermometers at 
intake pipes, or extrapolation methods. 

Here, we want to update the feedback law depending 
on intermediate disturbances w and newly obtained state 
estimates x. Hence, while shifting and iterating the prob- 
lem remains unchanged, a modified control JI N (x,w) is 
implemented instead of Pn{x, w) and we assume pZ N (x, w) 
to be instantly computable. 



For simplicity of exposition, we assume that a minimizcr 
u*(-,x,w) of (2) exists for each x e X and N e N. 
Particularly, this includes the assumption that a feasible 
solution exists for each x £ X and each w € W. For 
methods on avoiding the feasibility assumption in x we 
refer to Primbs and Nevistic [2000] or Grime and Pannck 
[2011]. Using the existence of a minimizer u*(-,x,w) G 
U x (x, w), we obtain the following equality for the optimal 
value function defined on a finite horizon 

Vn{x,w) = inf Jn(x,u,w). (3) 

In order to compute a performance or suboptimality index 
of the updated MPC feedback Ji N (x,w), we denote the 
closed loop trajectory by 

x(n + 1) = f(x(n),JI N (x(n),w(n)),w(n)) 
which gives rise to the closed loop costs 

oo 

V£»(x,w) :=^t(x(n),]i N (x(n),W(n)),w(n)). (4) 

Regarding the disturbances, we assume the following: 

Assumption 2.2. The disturbance is bounded by Aw > 
and the maximal state estimate deviation is bounded by 
Ax > 0, i.e. \\w\\ w < Aw and < Ax. Additionally, 

the nominal disturbance change is bounded by Aw, that 
is \\w(n+ l)|b ( „) < Aw. 

Due to the presence of unknown disturbances, asymp- 
totic stability cannot be expected. Therefore, we consider 
the concept of P-practical asymptotic stability. Similar 
to input-to-state stability (ISS), the convergence of the 
closed loop solution is characterized by comparison func- 
tions. Here, we call a continuous function p : R>o — > ]R>o 
a class /Coo-function if it satisfies p(0) = 0, is strictly in- 
creasing and unbounded. A continuous function /? : K> x 
M>o — > K>o is said to be of class ICC if for each r > 
the limit lim^oo (3(r,t) — holds and for each t > the 
condition (3(-,t) € /Ceo is satisfied. 

Definition 2. 3. Let A C X be a forward invariant set with 
respect to all possible disturbances satisfying Assumption 
2.2 and let P C A. Then x* e P is P-practically 
asymptotically stable on A if there exists (3 e JC£ such 
that 

\\x(n)\\ x *<P(\\x \\ x .,n)^ (5) 
holds for all xo € A and all n € No with x(n) £ P. 

Note that the ISS property can be shown for P-practical 
asymptotically stable systems by a suitable choice of the 
comparison function 7 € /C, cf. Chapter 8.5 in Grime and 
Pannck [2011]. Typically, the ISS property is shown via 
ISS Lyapunov functions, see, e.g, Jiang and Wang [2001] 
or Magni and Scattolini [2007]. Different from that, P- 
practical asymptotic stability can be concluded if there 
exists a suitable "truncated" Lyapunov function as shown 
in Theorem 2.20 in Grime and Pannck [2011]. 

Theorem 2.4- Suppose A C X is a forward invariant 
set with respect to all possible disturbances satisfying 
Assumption 2.2, P C A and x* e P. If there exist JC- 
functions ot\ , ct2 , 013 and a Lyapunov function V on S = A\ 
P satisfying 

aiflNU* + \\w\\w*) < V(x,w) < a 2 (|N|x* + MU*) 
V(x,w) > V(f(x,-p N (x,w),w),w) 
-a 3 (\\x\\x*) 



then x* is P-practically asymptotically stable on A. 
3. STABILITY 

In the literature, one usually uses terminal constraints 
or Lyapunov type terminal costs to guarantee the ISS 
property of the disturbed closed loop, cf., e.g., Diehl et al. 
[2005], Bock et al. [2007], Zavala and Biegler [2009]. Here, 
we consider the plain MPC formulation without these 
modifications. Instead, we suppose a relaxed Lyapunov 
condition to hold for the nominal case which additionally 
reveals a performance index of the nominal closed loop, 
cf. Lincoln and Rantzcr [2006] and Grime and Rantzer 
[2008]. Note that this condition is always satisfied if 
N is sufficiently large, see Alamir and Bornard [1995], 
Jadbabaie and Hauser [2005] or Grimm et al. [2005]. 

In order to prove a similar result in the disturbed case, 
one could modify the stage cost I to be positive definite 
with respect to a robustly stabilizable forward invariant 
neighbourhood of x* . Since the computation of this neigh- 
bourhood may be impossible, we choose the stage cost 
I to be positive definite with respect to x* only, that is 
ignoring the effects of disturbances on stabilizability. As I 
is typically much smaller close to x* than far away from 
the desired steady state, we may still expect the closed 
loop to converge to a neighbourhood of x*, i.e. P-practical 
stability of the closed loop. Similar to results shown in 
Grime and Rantzer [2008] we additionally obtain a bound 
for the degree of suboptimality. 

In order to show such a performance results, we assume 
the following: 

Assumption 3.1. There exists a set A C X containing x* 
such that 

\£(x,~p N (x,w),w) -£(x,hn(x, w),w)\ 

< L e (\\x\\ x + \\w\\ w ) 
\J N (x,Ji N (x,w),w) - J N (x,fi N (x,w),w)\ 

<Lj(\\x\\ x + \\w\\ w ) 

hold with Lipschitz constants Li and Lj for all tupels 
(x,x,w,w) with x, x £ A, w, w £ W satisfying Assump- 
tion 2.2. 

Under these conditions, we can prove the following: 
Theorem 3.2. Suppose Assumptions 2.2, 3.1 to hold and 
a given feedback [in ■ X x W — > U and a nonnegative 
function Vn ■ X x W — » to satisfy the relaxed 
Lyapunov inequality 

V N (x,w) > V N (f(x,(i N (x,w),w),w) + a£(x,fi N (x,w),w) 

(6) 

for some a £ (0, 1) and all x £ A and all w £ W. Suppose 
e > Lt(Ax + Aw) + (2LjAx + 3LjAw)/a and let L denote 
the minimal set which contains x* , is forward invariant 
with respect to all possible disturbances satisfying As- 
sumption 2.2 and Vn(x,w) > VN(f(x,[iN(x,w),w) + ae 
holds for all x £ A \ C and all w £ W. Furthermore define 
the modified costs 

-.. . (m&x{£(x, u, w) — e, 01 x £ A \ C ,„ N 

*(*,«,«,):= ( xe£ " (7) 

and a := mf{VN(f(x,fj,N(x,w),w),w) | x G A\C,w € 
W}. Then for the modified closed loop cost 



V£> N {x,w) := £{x{n), jjL N {x{n),w(n)),w{ri)) 

n=0 

we have 

a V^ N (x, w) < Vn (x, w) - a <Vca (x, w) - a (8) 
for all x £ A. 

Proof: Consider xo £ A. Let uq £ No be minimal with 
x(no + 1) £ C and set n := oo if this case does not occur. 
Reformulating (6) we obtain 

a£(x,fi N (x, w),w) < V N (x,w) - V N {f{x,n N {x,w),w),w). 

Now we incorporate the effects of disturbances and sensi- 
tivity updates of the control using Assumption 3.1 which 
gives us 

a£(x(n),J2 N (x(n),w(n)),w(n)) 

< a£(x(n), fiN(x(n), w(n)), w(n)) 

+ aL e (\\x(n)\\ x{n} + \\w(n)\\ w(n) ) 

< V N (x(n),w(n)) + aLi (||s(n)|| x( „) + \\w(n)\\ w{n) ) 

- V N (f(x(n),fj, N (x(n),w(n)), w(n)),w(n)) 

< V N (x(n),w(n)) + aLi (\\x(n)\\ x(n) + \\w(n)\\ w{n) ) 

- V N (f(x(n),fi N (x(n),w(n)),w(n)),w(n + 1)) 
+ Lj\\w(n + l)\\ w ( n ) 

< V N (x(n),w(n)) - V N (x(n+l),w(n+l)) 

+ aL e (\\x(n)\\ x(n) + \\w(n)\\ w(n) ) 
+ Lj\\w(n+ l)|U (n ) 

< V N (x(n),w{n)) - V N (x(n + l),w(n+l)) 

+ (aL e + Lj) (\\x(n)\\ x(n) + \\w(n)\\ w{n} ) 
+ Lj (\\x(n + l)|U ( „+i) + \\w(n + l)|U(„+i)) 
+ Lj\\w(n+ l)\\ w ( n) 
Hence, using boundedness from Assumption 2.2 reveals 
a£(x(n) , ~p N (x(n),w(n)),w(n)) 
< V N (x(n),w(n)) - V N (x{n + l),w{n + l)) 
+ (aL e + 2Lj)Ax + (aL e + 3Lj)Aw. 

Since we have that Vn (x, w) > Vn(/(x, [in(%,w),w),w) + 
(aL t + 2Lj)Ax + (aL e + 3Lj)Aw holds for all x £ A \ 
C, we obtain Vn(x,w) > Vn(/(x,]I n (x,w),W),w) for all 
x £ A \ C and all w satisfying Assumption 2.2. Using this 
fact, the bound on e and the definition of I in (7) we have 

a£(x(n) , JI N (x(n) , w(n)), w(n) ) 
= max {a£(x(n), ~p N (x(n), w(n)), w(n)) — ae, 0} 
< V N (x(n),w(n)) - V N (x(n + l),w(n+l)). 
For n > uq + 1 the invariance of C gives us x(n) £ C and 
hence £(x(n),JI N (x(n),w(n)),w(n)) = 0. Additionally, 
since a is the minimal cost after entry in C, we have 
Vn(x(u), w(n)) > g for all n < n + 1. Now we can sum 
the stage costs over n and obtain 

K 

ot^£(x(n),Ji N (x(n),w(n)),w(n)) 

n=0 

K _ 

= a^£(x(n),-p N (x(n),w(n)),w(n)) 

< V N (x(0),w(0)) - V N (x(K + l)w(K + 1)) 

< V(x(0),w(0)) - a. 



where Kq := min{K, no}. Using that KeN was arbitrary 
we can conclude that (Vn{x(0),w(0)) — a)/a is an upper 

bound for V£ N (x,w). □ 

The definition of A in Theorem 3.2 is implicit, yet an ap- 
proximation of A can be obtained via techniques presented 
in Grimm ct al. [2005] or Grime and Rantzer [2008]. 

In addition to the previous performance estimate, im- 
practical asymptotic stability can be shown as follows: 

Theorem 3.3. Suppose the conditions of Theorem 3.2 hold 
and additionally there exist JC functions ai, a 2 such that 

ai(||z|U* + IMU*) < V(x,w) < a 2 (\\x\\x* + \\w\\ w *) 
holds for all x € A \ P with P = C. Then x* is P- 
practically asymptotically stable on A. 

Proof: Follow directly from the definition of ~£ in Theorem 
3.2 and the property V(x, w) > VN(f(x,JI N (x,w),w),W) 
shwon in the proof of Theorem 3.2. □ 

Note that the stability result holds for any update on the 
feedback law ~p N and is not specific to how this update is 
obtained. 

Remark 3.4- Condition (6) can be relaxed to 

V N (x,w) - V N (f(x,fj, N (x, w),w), w) 
> mm{a£(x, Hn(x, w),w) — e, £(x, Hn(x, w), w) — e} 

to additionally allow the nominal system to be P- 
practically stable only, cf. Griine and Rantzer [2008]. In 
this case, we have to modifiy the lower bound on e to 
((Li + 2Lj)(Aw + Ax) + e)/a and require Vn(x,w) > 
Vjv(/(x, Hn{x, w), w) + 2Lj(Aw + Ax) + e to hold for all 

x e A \ C. 

After stating stability for an abstract update law of ~p, N , 
our goal in the next section is to verify the required 
Assumption 3.1 for a particular updating strategy. 

4. SENSITIVITY THEORY 

Sensitivity analysis has been analysed extensively for the 
case of open loop optimal controls, see, e.g, Grotschel et al. 
[2001], and has become a rather popular control method. 
In the MPC feedback context, Zavala and Biegler [2009] 
analysed the impact of sensitivity updates on stability 
of the closed loop in the presence of stabilizing terminal 
constraints and Lyapunov type terminal costs. Here, we 
use an identical update 



' du* 



u(-,x,w) := u*(-,x, w) + 



{■,x,w) 
-(;x,w) 



x(-)-x(-) 
w(-) - w(-) 



(9) 

and set ]J N (x,w) — u(0,x,w) but consider the plain MPC 
case described in Section 2. Since the solution of the 
underlying optimal control problem and the computation 
of the sensitivities require some computing time, such an 
approach is typically implemented in an advanced step 
setting, see also Findeisen and Allgower [2004]. The idea 
of the advanced step is to precompute an open loop 
control u*(-,x,w) and the sensitivities for a 

future time instant while the current sampling period 
evolves. Then, once the time instant is reached at which 



the control is supposed to be implemented, newly obtained 
measurements of x and w are used to update the control 
according to (9). An algorithmic implementation of this 
idea is shown below. 

MPC Algorithm 

(1) Obtain measurement of x(n) and w(n) 

(2) Update control u(-,x(n),w(n) (e.g. via (9)) and apply 
control JI N (x(n),w(n)) = u(0,x(n),w(n) 

(3) Predict x(n+l) using dynamic (1) together with x(n), 
Ji N (x(n),w(n)) and w(n) 

(4) Compute u*(-,x(n + l),w(n)), set n := n+l and goto 
step (1) 



The mathematical foundations of such a control update 
are given in Fiacco [1983]. For notational convenience, we 
adapted these results to the considered MPC case. 

Theorem 4-1- Suppose that / and £ are twice continuously 
diffcrentiablc in a neighbourhood of the nominal solution 
u*(-, x, w). If the linear independence constraint qualifica- 
tion (LICQ), the sufficient second order optimality con- 
ditions (SSOC) and the strict complementarity condition 
(SCC) are satisfied in this neighbourhood, then we have 
that 

• u*(-,x,w) is an isolated local minimizer and the 
respective Lagrange multipliers are unique, 

• for (x, w) in a neighbourhood of (x, w) there exists 
a unique local minimizer u*(-,x, w) which satisfies 
LICQ, SSOC and SCC and is differentiable with 
respect to x and w, 

• there exist a Lipschitz constant Lj such that 

| J N (x,~p N (x, w),W) - J N (x,^ N (x 7 w),w)\ 

<Lj(\\x\\ x + \\w\\ w ) (10) 

holds and 

• there exist a Lipschitz constant L u such that for 
the updated control u(-,x,w) from (9) the following 
estimate holds: 

d u (u(x,w),u*(x,w))) < L u (\\x\\ x + \\w\\ w ) 

Using this result, Assumption 3.1 can be verified. 

Proposition 4-2. If there exists a set A C X containing 
x* such that conditions of Theorem 4.1 hold for all 
tupcls (x, x,w, w) with x,x € A, w,w € W satisfying 
Assumption 2.2, then Assumption 3.1 holds. 

Proof: Follows directly from (10) and the fact that dif- 
ferentiability of £ implies the existence of a local Lipschitz 
constant Lg. □ 

Note that in order to apply both Proposition 4.2 and 
Theorem 3.2 the set A is not necessarily large, a fact 
that otherwise may exclude such an approach in the 
presence of state and control constraints. In particular, 
Theorem 4.1 requires the open-loop control structure to 
remain unchanged in a neighbourhood of the optimal 
solution despite changes in x and w. Most importantly, 
this property is only required locally. On a larger scale, the 
closed-loop control structure may change since we allow 
for an intermediate reoptimization. Hence, despite the fact 
that the update formula (9) is restricted to a certain 
neighbourhood of the open-loop solution, the MPC update 



approach using sensitivities is only locally restricted to 
that particular neighbourhood, i.e. for each visited closed- 
loop state this neighbourhood and the respective control 
structure may change. 

5. NUMERICAL RESULTS 

To illustrate our results, we consider a halfcar model given 
from Speckert et al. [2009], Popp and Schiehlen [2010] with 
proactive dampers given by the second order dynamics 

mixi = mig + / 3 - h 
m 2 x 2 = m 2 g 4-/4-/2 
7713X3 = m 3 g - / 3 - f 4 
Ix 4 = cos(x 4 )(bf 3 - af 4 ) 
where the control enters the forces 
fi = k i (x l - Wi) + di(±i - Wi), i = l,2 
/3 = ^3(^3 — xi — bsm(x 4 )) + ui(x3 — ±1 — bx 4 cos(a;4)) 
f 4 = k 4 (x 3 — x 2 + asin^)) + u 2 (x 3 — x 2 + ax 4 cos(x4)) 
see Fig. 1 for a schematical sketch. Here, x\ and x 2 denote 




w(t — A) w(t) 



Figure 1. Schematical sketch of a halfcar subject to road 
excitation to 

the centers of gravity of the wheels, 23 the respective 
center of the chassis and x 4 the pitch angle of the car. 
The disturbances Wi, w 2 are connected via wi(t) = w(t), 
w 2 (t) = w(t — A). The remaining constants of the halfcar 
are displayed in Tab. 1. For this problem, we apply MPC 



name 


symbol 


quantity 


unit 


distance to joint 


a, b 


1 


777 


mass wheel 


mi , 7772 


15 


kg 


mass chassis 


7«3 


750 


kg 


inertia 


/ 


500 


kg m 2 


spring constant wheels 




2 ■ 10 5 


kN/m 


damper constant wheels 


di,d,2 


2 ■ 10 2 


kNs/m 


spring constant chassis 


k3,k 4 


1 • 10 5 


kN/m 


gravitational constant 


9 


9.81 


m/s 2 



Table 1. Parameters for the halfcar example 



cost functional 

JV-l 

J N (x,u,w) := ^2 VRFn(k) + [i A F A (k) 
k=o 

following ISO 2631 with horizon length N = 5. The 
handling objective is implemented via 



F R{k):= J ^ [fci (si (g - m (*)) + di (M (t) - vn m -F ± y 

kT 

+ / [fcg (a-a (t) - t»2 (t ) ) + d 2 ( ±2 (t) - w 2 (t) )] - F 2 ^ 2 ^ 

with nominal forces F\ — (a - g ■ (mi + 777,2 + m 3,))l ( a + b), 
F 2 = (b ■ g ■ (mi + fn-2, + m 3))/( c ' + b) whereas minimizing 
the chassis jerk 

(k+l)T 

F A (k) := J (m 3 x 3 (t)) 2 dt 

kT 

is used to treat the comfort objective. Both integrals are 
evaluated using a constant sampling rate of T = 0.1s 
during which the control are held constant, i.e. the control 
is implmcntcd in a zero-order hold manner. Additionally, 
the control constraints U = [0. 2kNs/m, 5kNs/m] 2 limit 
the range of the controllable dampers. Within the MPC 
scheme outlined in Algorithm 4, we compute the nominal 
disturbance ttj(-) and the corresponding derivates from 
road profile measurements taken at a sampling rate of 
0.002s via a fast Fourier transformation (FFT). Both the 
states of the system and the road profile measurements 
are modified using a disturbance which is uniformly dis- 
tributed in the interval [-0.025,0.025]. 

Jerk 
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Figure 2. Chassis jerk for MPC with (blue) and without 
sensitivity update (red) 

As expected, the updated control law shows a better 
performance which can not only be observed from Fig. 
2, but also in terms of the closed loop costs: For the 
chosen industrial road data we obtain an improvement 
of approximately 8.2% using the sensitivity update (9). 
Although this seems to be a fairly small improvement, 
the best possible result obtained by a full reoptimization 
reveals a reduction of approximately 10.5% of the closed 
loop costs. 

We like to note that due to the presence of constraints it is 
a priori whether the conditions of Theorem 4.1 hold at each 
visited point along the closed loop. Such an occurrance can 
be detected online by checking for violations of constraints 
or changes in the control structure. Yet, due to the 
structure of the MPC algorithm, such an event has to be 
treated if one of the constraints is violated at open loop 



time instant k = 1 only which was not the case for our 
example. 

6. CONCLUSIONS AND OUTLOOK 

We considered MPC without stabilizing terminal con- 
straints or Lyapunov type terminal costs in the presence 
of disturbances. For this setting, we presented stability 
and performance results for abstract control updates. The 
required assumptions are shown to be fulfilled if sensitivity 
updates on the control can be performed. 

Future research concerns the analysis of other types of con- 
trol updates such as realtime iterations which have been 
considered for MPC with stabilizing terminal constraints 
and costs. Using the presented result, we hope to obtain a 
unified stability and performance analysis of MPC in the 
presence of disturbances and respective control updates. 
Regarding the halfcar example, the influence of other 
types of interpolation/extrapolation methods on the MPC 
performance such as road profiles based on statistical data 
will be analysed. 
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